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Abstract— A question that often comes up in the context 
of an IRA is “how is the transient energy from the pulser 
radiated in space”? Of course the electromagnetic fields 
(both E and H), the power density and the energy density 
have their maximum on the bore sight. Power pattern is a 
well-defined frequency domain concept, but it is a 
cumbersome descriptor for hyperband antennas such as an 
IRA because of the multitude of frequencies involved. In 
this paper we explore the concept of an energy pattern 
which holds good both in time and frequency domains. An 
energy pattern is useful in visualizing where the transient 
energy provided to the IRA is going. It is further noted that 
the energy and power patterns are identical for a CW 
antenna, while they can be vastly different for pulsed 
antennas. 


Index Terms—Impulse Radiating Antenna, Pulsed 
Antenna, Radiation Characteristics, Power Pattern, 
Energy Pattern 


I. INTRODUCTION 

An easy way to remember the performance parameters of 
an antenna is through an acronym “BRIDGE”. The 
individual letters denote: 
Beam width 
Radiation pattern 
Input impedance 
Directivity 
Gain 
Effective area 


All of these parameters are well defined in the frequency 
domain, and are functions of frequency. In the context of 
pulsed antennas where many frequencies are simultaneously 
fed into the antenna, the use of frequency dependent 
parameters is useful, but cumbersome. 


Let us consider a reflector type of an IRA as an example 
of a hyper band [1] antenna [2, 3, and 4]. The radiation 
pattern of an IRA is a strong function of frequency as 
reported in [5]. The lower frequencies of the input pulse have 


lower gain and large beam widths, while the higher 
frequencies have a higher gain and smaller beam widths. In 
the next section, we consider an energy pattern of the IRA 
which is a simple and unique descriptor. It is indicative of 
how the input energy is radiated into all of space. 


Il. ENERGY PATTERN 


Let us denote the energy pattern of the IRA by U(0,@) in 


the far field or Fraunhofer zone. This quantity is measured in 
Joules/steradians. Figure 1 shows the geometry and Cartesian 
and spherical sets of coordinates with origin at the center of 
the radiating aperture. The diameter and the focal length of 
the reflector are denoted by D and F. The Fraunhofer zone is 
known to begin at an axial distance given by [6] 
D ? V peak 


fp En ; = > 1 
far 2CT p, Tr (=) ( ) 
dt peak 


where c is the speed of light, D = diameter of the reflector, 
and T,,, is the maximum rate of rise of the voltage wave 
launched onto the reflector, which is not necessarily the rate 
of rise of the transient source waveform. There could be 


some risetime degradation in transporting the voltage pulse 
from the source to the wave launch. 


y 
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Figure 1. The IRA geometry and coordinate system. 


Let Er (R,0,@,t) and E(R,0,0,0) denote the far 


electric field at an arbitrary location in time and frequency 
domains respectively. The energy pattern can be defined as 
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Joules/Steradian (2) 


Note that equation (2) represents an energy density and not 
the total radiated energy. In equation (2), Zọ is the 
characteristic impedance of free space and the value of r 
needs to be in the far field satisfying equation (1). It is noted 
that energy content of a signal is the same whether one 
computes it in time domain or frequency domain, as per 
Parseval’s Theorem [7]. It is also observed while the 
electromagnetic fields depend on the radial distance r, but the 
energy pattern in the far field is independent of r, simply 


because the electric field is falling off like r. We can re- 
write equation (2), in short hand notation as 
U -fr [Eo] it (Joules /steradian) (3) 
Z 


0 o 


In equation (3), E(t ) denotes the total electric field or the 
vector sum of all non-zero components. At an arbitrary point 
(r,0,@) in the far field, the electric field has a vanishing 
radial field, and there are only @ and ø components. The 


usual relationships between the rectangular and spherical 
components are applicable here, if one chooses to compute 
the Cartesian complements of electric field. 


E.=0 
E, = E ,cos(@) cos(g) + E, cos(@) sin(g)—E, sin (0) 
E, =—£, sin(g) + E,cos(9) (4) 


E, = E, cos (0) cos (p) — E, sin (p) 
E, = E, cos (8)sin (9) + E cos (o) 
E, =- Esin (0) 


The square of the magnitude total electric field in time 
domain is given by 


REXA =[E} + E] =[E + E+E] © 


Of course, in practice, the integration in equation (3) is 
carried out over time limited or band limited ranges. 


HI. ESTIMATION OF FAR FIELD BY APERTURE 
INTEGRATION METHOD 


The electric field at any arbitrary observation point in the 
far field can be estimated by the method of aperture 
integration that takes into account the presence of the feed 
arms. Initially the field on the aperture plane is computed. 
The radiated field can then be computed by an integration of 
the magnetic current (tangential electric field) on the 
aperture. The relevant equations are given below [8], with 
reference to Figure 1. 


The magnetic current on the aperture is given by 


M=-2 Í xĒ (6) 


z aperture 


and the electric field components in terms of the magnetic 
current on the aperture are given by 











1 7 1+ jkr y oe 
E, = rps ffe z) M, 3 e™ dx'dy 
1 r 1+ jkr — jkr ' $ 
E, =; ffe z')M, =e H dx" dy 
1 BS ee 
E,=7—J[to-y)M,-(@-x)M,I—— &™ dr dy 
(7) 


From the estimation of all components of the electric field 
in frequency domain, the total transient electric field is found 
by Fourier inversion and then used in equation (3) to get the 
energy pattern. 


Some comments about the aperture integration method 
outlined above are in order. 


The aperture integration method for a hole in an infinite 
screen is quite accurate in the direction of the main beam, but 
the accuracy deteriorates as the angle @ increases. This is 
due to the various approximations used to obtain an 
analytical solution to this problem. In the EM model used 
here for the dish, the tangential E and H-fields are normally 
needed over the aperture, and are assumed to be zero 
elsewhere on the aperture plane away from the dish. We use 
only the tangential E-field distribution, however, and assume 
the corresponding TEM value for the tangential H-field (and 
this results in an equivalent magnetic current of M = 2 E, 
over the aperture). 

For this idealized aperture source, the radiated field in the 
E-plane when @= (plus or minus) 90 degrees is zero due to 
geometrical projection reasons. In the H plane, the radiated 
field does not vanish for @ = (plus or minus) 90 degrees, 
because the equivalent source is radiating in free-space and 
this observation point is broadside to the magnetic current 
source. If there were an infinite, perfectly conducting plate in 
the plane of the aperture, there the field in the H-plane would 
go to zero at + - 90 degrees due to the boundary condition of 
Eran = 0 on the screen. We also observe that the aperture 
integration model will not give good results on the back-side 
of the aperture plane. In fact, the radiation pattern from the 
equivalent sources is the same in the back as in the front. 


In summary, the aperture integration model gives good 
responses in the forward direction, but not in the backward 
direction. The numerical errors also increase as the 
observation point approaches the plane of the aperture. 


IV. ENERGY PATTERN OF THE PROTOTYPE IRA 


In this section, we consider the prototype IRA and 
evaluate its energy pattern. A photograph of this prototype 
IRA is shown in Figure 2. 
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Figure 2. Photograph of the Prototype IRA. 


The parameters of the prototype IRA are summarized in 
Table 1, where we have separately listed the geometry of the 
reflector and the pulser parameters. Having outlined the 
salient features of the prototype IRA, we now proceed to 
compute the radiated fields and the energy pasterns. 

The excitation voltage of the IRA is the total plate voltage 
situated at the focal point of the dish, and is denoted as Vo (^). 
As outlined in [6], one possible representation of a fast-pulse 
is given by the expression 


Vo(0 = Voeu: (I+ p ose =) u(-( -1))e| teas we =) u(t -) 


(8) 
TABLE 1. Antenna and pulser parameters of the 
prototype IRA 


(a)Antenna Details 
































Parameter Value 
Reflector diameter D 3.66.m 
Focal length F 1.22m 
F/D 0.33 
Number of arms 4 
Arm configuration 90 deg 
Impedance Z, 200 Ohms 
Geometrical factor 0.53 
T, eo Zin / Zo 
Polarization Linear/Vertical 











(b) Pulser Details 





Parameter 


Value 





Peak voltage Vp 


~+60kV ~ 120 kV 





Peak rate of rise 

















(£) ~ 1.2 x 10" V/s 
dt JP 
Maximum rate of rise 
pus Vp ~ 100 ps 
mr dV 
dt), 
Pulse repetition frequency ~200 Hz 
PRF 
Pulse decay time t} ~ 20 ns 
5 
Duty cycle T = t} PRF ~4X10 
y? ~72 MW 
Peak power P, = —> 
Z 


in 








Average power 





~72MW x4x 10° 
~ 288 Watts 
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In equation (8), the term erfc(.) denotes the 
complementary error function and u(.) is the unit step 
(Heaviside) function. For the calculations here, the following 
parameters are used: 
© V eak = 120 (kV) (peak value of the transient 
waveform) 

e T =0.006 (amplitude adjustment factor) 

e £= 0.005 (fall-time coefficient) = rise time / fall 
time = 100ps / 20 ns 


e T, =100 (ps) (waveform rise time) 


e T, = 1.0 (ns) (time shift) 


Figure 3 presents the waveform V,(t) for late and early 


times, respectively, 
derivative dV/,(t)/dt 


magnitude of the waveform is provided in Figure 5. 


and Figure 4 illustrates the time 
of the waveform. The spectral 
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Figure 3. Plot of the excitation waveform V(t). 
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Figure 4. Plot of the time rate of change of the excitation 
waveform dV/,(t)/dt. 
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Figure 5. Plot of the spectral magnitude of the excitation 
waveform. 


The calculation of the radiated far-field response of this 
antenna is done using an aperture integration method, which 
is described in [9]. The TEM field illuminating the dish is 
approximated and this forms the aperture field over the dish 
which is integrated to yield the impulse portion of the 
response. In addition, the radiation from the feed line 
currents and the source itself is determined and this provides 
the pre-pulse contribution to the field. 


Overlay plots of the principal components of the transient 
E-fields for this antenna are presented in Figure 6 for the H- 
plane and in 7 for the E-plane. 


As a check of these results, Giri [6] provides the following 
estimation of the peak value of the impulse portion of the 
normalized radiated E-field on axis (0, ø) = 0°, 0° : 


D42 V 
sP (=) V2 [eee = evan 
peak 


TE seak 
4rcf,\ dt ARCH T mr 





(9) 


In this expression, OV,/Ot|_ =~ 1.184*10!5 (V/sec) 





peak 
from Figure 4 and f, = 1.06. This provides an estimate of 
rEpeak = 1.534 MV, which is shown in the figures and 
agrees well with the present aperture integration calculations. 
In these responses, it should be noted that the large 
impulse-like response occurs only in the boresight direction. 
It occurs in a very narrow beam around (8, ø) = 0°, 0°, and 
even at a value of 0 = 1° there is a substantial reduction of 
the peak. Moreover, the waveforms become more dispersed 
in time as the angle @ increases. It is also noted that the 
waveforms in the H-plane and in the E-plane are different. 
Although the E£, and E,, field components are not plotted 


here, it suffices to note that they are both either zero or 
negligible compared with the principal Ey component. 
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Figure 6. Overlay of the principal (Ey) field components 


for various values of the angle 8, computed in the horizontal 
plane (the H-plane) in the far zone of the IRA. (for 
120k V/100ps/20ns pulse excitation). 


The spectral responses of the Ey field components in the 


H-plane and E-plane are also of interest. These are plotted as 
overlay plots for different values of the observation angle 0 
in Figure 8 and Figure 9 respectively. 


For the driving waveform of Figure 3 into an impedance 
of 200 Ohms (the net IRA impedance), the total energy 
delivered to the antenna can be estimated to be 


1% 


N-1 
U, = f v2 Odi s ALS (V2) =0.727 Joules (10) 
2004 200 — 


Of interest is the spatial distribution density (in 
Joules/steradian) of the radiated energy from the IRA. 
Equation (3) is evaluated in both of the observation planes by 
and the results are presented in Figure 10 and Figure 
11Figure for different scales for the y-axis. 
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Figure 7. Overlay of the principal (Ey) field components 


for various values of the angle @, computed in the vertical 
plane (the E-plane) in the far zone of the IRA. (for 
120kV/100ps/20ns pulse excitation). 
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Figure 8. Plot of the spectral magnitudes of the principal 
(Ey) field components computed in the horizontal plane (the 


H-plane) in the far zone of the IRA. (120kV/100ps/20ns 
pulse excitation.) 
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Normalized E-field Spectrum rEy (kV/Hz) 
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Figure 9. Plot of the spectral magnitudes of the principal 
(Ey) field components computed in the vertical plane (the E- 


plane) in the far zone of the IRA. (120kV/100ps/20ns pulse 
excitation. ) 


In this calculation, all vector components of the calculated 
E-field are used, even though x and y-components are very 
small. It is clear that the IRA sends more energy in the 
boresight direction. It is important to note that this radiated 
energy pattern will be different for other excitation sources, 
because it depends on the frequency content of the excitation 
voltage. 
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Figure 10.Plot of the radiated energy pattern in the 
horizontal (H) and vertical (E) planes for the IRA with the 
voltage excitation of Figure 3. (Log scale on the ordinate.) 
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Figure 11.Plot of the radiated energy pattern in the 
horizontal (H) and vertical (E) planes for the IRA with the 
voltage excitation of Figure 3. (Linear scale on the ordinate.) 


V. CALCULATIONS FOR THE PROTOTYPE IRA WITH A 
DAMPED SINE WAVE EXCITATION 


Because the radiated energy pattern depends on the 
waveform exciting the IRA, it is interesting to consider an 
alternate excitation function. In this section we examine a 
damped sine waveform with a center frequency of fọ = 1 


GHz and 
a= 2 x 108 (1/sec). The peak amplitude of this waveform is 


adjusted to provide the same value of the previous pulse 
waveform, namely Vmax = 120 kV. 


exponential damping constant 


Figure 12 illustrates this waveform, and its derivative is 
shown in Figure 13.. The spectral magnitude is provided in 
Figure 14. For this waveform, the resulting transient 
responses in the H- and E-planes are shown in Figure 15 and 
Figure 16, and the spectral responses are in Figure 17 and 
Figure 18.and Figure. It is clear that the responses for this 
excitation are considerably different from those of the fast 
pulse. 

Note that the estimated peak value of rEpeak = 1.009 MV 
from equation (9) agrees well with the computed results. It is 
also worth noting that the low amplitude early time ringing 
in the waveforms of Figure 15 and Figure 16 is not a FFT 
problem, but rather, it is the pre-pulse response from the 
IRA. 


Vo (kV) 


Figure13. Plot of the time rate of change of the excitation 
waveform dV,,(1)/dt. 
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Figure 14.Plot of the spectral magnitude of the excitation 
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Figure 12.Plot of the excitation waveform V/(f). 
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Figure15. Illustration of the principal (Ey) field component 


computed in the horizontal plane (the H-plane) in the far 
zone of the IRA. (120kV/ 1 GHz Damped sinusoidal pulse 


excitation.) 
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Figure16 Illus 


computed in the 
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tration of the principal (Ey) field component 


vertical plane (the E-plane) in the far zone 


of the IRA. (120 kV/ 1 GHz Damped sinusoidal pulse 


excitation.) 
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Normalized E-field Spectrum rEy (kV/Hz) 
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Figure17. Plot of the spectral magnitudes of the principal 
(Ey) field component computed in the horizontal plane (the 


H-plane) in the far zone of the IRA. (120 kV/ 1 GHz 
Damped sinusoidal pulse excitation.) 
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Figurel8.. Plot of the spectral magnitudes of the 
principal (Ey) field component computed in the vertical 


plane (the E-plane) in the far zone of the IRA. (120 kV/ 1 
GHz Damped sinusoidal pulse excitation.) 


Figures 19 and 20 present the radiated energy pattern for 
this damped sine excitation. It is noted that there is more 
energy radiated in the main beam for this excitation than for 
the fast pulse. In addition, there are significant side lobe 
variations in the pattern for off-axis directions. This was not 
the case for the fast-pulse excitation, because this former 
excitation contains many frequencies with similar amplitudes 
and the side lobes tend to wash out in the transient response. 


This is not the case for the damped sine wave, which has a 
rather narrow range of significant frequencies in the 
spectrum. 


VI. SUMMARY 


We have explored the concept of energy patterns 
(measured in Joules/steradian) for a pulsed antenna such as 
the IRA. In the case of CW antennas, the radiated power and 
energy patterns are the same. This is not the case for pulsed 
antennas. The radiated power pattern for a pulsed antenna is 
a strong function of frequency and can be computed for 
various frequencies. 


In this paper, we have considered the prototype IRA and 
estimated its energy pattern for two different inputs with the 
same peak voltage amplitude of 120 kV. One input is a fast 
rising (100ps) - slowly decaying (20ns) mono-polar pulse, 
while the second is a 1 GHz damped sinusoidal voltage, 
which is bipolar. 


The transient energies of the two input voltages have 
vastly different frequency components. The fast pulse has 
frequencies extending from DC to a few GHz, while the 
damped sinusoidal input is a moderate band source centered 
at 1 GHz. As a consequence of this, the energy patterns of 
the same prototype IRA are considerably different for the 
two input voltages. 


In summary, the energy patterns are useful in visualizing 
where the transient energy provided to the IRA is being 
radiated in front of the antenna. 


This paper is an adaptation of [10]. 
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